Short-time scaling in the critical dynamics of an antiferromagnetic Ising system with 

conserved magnetisation 
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We study by Monte Carlo simulations the short-time exponent in an antiferromagnetic Ising 
system for which the magnetisation is conserved but the sublattice magnetisation (which is the order 
parameter in this case) is not. This system belongs to the dynamic class of model C. We use nearest 
neighbour Kawasaki dynamics so that the magnetisation is conserved locally. We find that in three 
dimensions 8 is independent of the conserved magnetisation. This is in agreement with the available 
theoretical studies, but in disagreement with previous simulation studies with global conservation 
algorithm. However, we agree with both these studies regarding the result 6c 7^ 9 a- We also find 
that in two dimensions, 9c = 9a- 



In equilibrium statistical physics, universal scaling laws 
are observed close to a critical point where the correla- 
tion length diverges. Dynamical systems also exhibit a 
universal scaling behaviour in the long time regime. Dy- 
namic universility classes are characterised by the dy- 
namic exponent, which connects the divergences in space 
and time. Typically, in a magnetic system, the finite size 
scaling form of a physical observable 0(t, r, L) is given 
by, 



0(t, r, L) = b- x O(b~ z t, b^ v r, b~ x L) 



(1) 



where t = (1 — T/T c ) is the deviation from the critical 
temperature T c , b the scaling factor, v and x the static 
critical exponents, z the dynamic exponent and L is the 
linear size of the system. 

Some time ago, it was found that there exists a uni- 
versal behaviour in short-time regime as well §. If a 
magnetic system is quenched from a high temperature to 
its critical temperature T c , with initial order parameter 
equal to mo, then universal scaling behaviour is observed 
in a macroscopic short-time regime : 

M(t, t, L, m ) = b~ p/v M{b- z t, b 1/u r, 6 -1 L, b x °m Q ) (2) 

where M is the order parameter, t is time, (3 is the criti- 
cal exponent associated with the order parameter and Xq 
is a new exponent associated with the short-time effect. 
xq is the scaling dimension of mo. At the critical point, 
for small values of mo, 



M(i, m ) ~ m t e 



(3) 



with 8 = (xq — /3/v)/z. is a new exponent, not related 
to any previously known static or dynamic exponent [jjj . 
The phenomenon is universal as does not depend on 
microscopic details, updating schemes etc. It is observed 
in a macroscopic short time regime. 

A positive value of 9 will indicate that the magneti- 
sation (order parameter in general) will first increase in 
time and later show the conventional power law decrease. 
It was observed in many systems that such an increase 
indeed occurs. The exact cause of this behaviour is not 
yet known very clearly. 



Systems with different critical dynamics have been 
classified as models A, B, C etc [Q. Model A has no con- 
servation, in model B the order parameter is conserved 
and in model C, the non-conserved order parameter is 
coupled to a non-ordering conserved field. These dy- 
namic classes are distinguished by the different values 
of the dynamic exponent z. 

Using field theoretical methods, it was shown in 
that short time behaviour exists in Model A. Extensive 
numerical studies for calculating in the different dy- 
namical classes have also been made in the recent years 
and accurate estimates of in model A are available in 
different dimensions ||. The results confirm the qual- 
itative behaviour predicted by the theoretical analysis. 
Systems belonging to the class of model B, do not show 
any short-time effect. However, numerical studies with 
a globally conserved order parameter in two dimensions 
indicate that there could be a short time effect Q] . 

In model C, field theoretic techniques have again shown 
that there is a universal short-time behaviour Q . Short 
time effect in a semi infinite system belonging to model C 
has also been studied . There has been some recent nu- 
merical studies of the short-time scaling in model C with 
global conservation . These numerical studies show the 
novel result that the short-time exponent depends on the 
value of the globally conserved magnetisation - a feature 
not obtained at least in the theoretical study with local 
conservation. However, this is not very surprising, as sys- 
tems with global conservation in general show different 
behaviour compared to the ones with local conservation 
|^| |l0| as far as the long time behaviour is concerned. In 
model C and model B for example, the behaviour with 
global conservation becomes model A-like. 

However, before any conclusive statement about the 
difference in short-time behaviour in model C with local 
and non-local conservation can be made, it is necessary 
to obtain the numerical estimates of the exponents where 
the coupling field is locally conserved. This is because the 
numerical estimates and field theoretic results could be 
quite different below the upper critical dimension. In this 
communication, we report the estimates of the short-time 
exponents in both two and three dimensions with several 
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values of the conserved density from simulation studies 
with local conservation. 

We have taken an antiferromagnetic Ising system 
where the magnetisation too is kept constant. This model 
has been studied numerically previously to estimate the 
dynamic exponent z JHJ and the short-time exponent 9 
with global conservation . The order parameter in this 
system is the staggered magnetisation m s which is not 
conserved. Initially the system is allowed to have short 
range correlations only (since we quench it from a high 
temperature to T c ) and the initial m s (t = 0) is kept at a 
low value. Hence the system is prepared with two con- 
straints: both mo and m s are kept fixed. This could be 
done by keeping a staggered field, but we simply keep it 
at the desired value by appropriately flipping spins in a 
random configuration (the method is analogous to what 
is called a sharply prepared state in ||). The antifer- 
romagnetic interactions are between nearest neighbours. 
The lattices are hypercubic. Periodic boundary condi- 
tion along one direction and helical boundary condition 
along the other directions are used. Helical boundary 
condition with antiferromagnetic nearest neighbour in- 
tercations demand that the linear size along (d — 1) di- 
rections should be odd where d is the spatial dimension. 

We use the Kawasaki dynamics where the opposite 
spins on nearest-neighbouring sites are exchanged, thus 
keeping mo a constant. The exchange between local 
neighbours ensure that the conservation is local. The 
spins are updated sequentially and a sweep through the 
entire system is equivalent to one unit of time (one Monte 
Carlo step). In both two and three dimensions, we fix too 
at several different values. Since the system is quenched 
to the critical temperature T c (too), we first estimate the 
value of T c (too) (Table 1) for a large system following the 
method used in . We then estimate 9 for each value of 
to in both two and three dimensions in smaller lattices. 
For this estimate, one should measure the initial slope of 
m s (t) vs t curve for different values of m s (0) and take 
the limit m s (0) — > 0. 

The results are as follows. In two dimensions, the 
short-time exponent 9 has been found to be 

9 = 0.18 ±0.01 (4) 

and this value has no measurable dependence on the ini- 
tial (small) value of sublattice magnetisation m s (0) (Fig. 
1) and on the conserved magnetisation to (Fig. 2). Also, 
in contrast with the simulation studies with global con- 
servation algorithm Q, it is the same for model A and 
C, upto the accuracy of the present study. In three di- 
mensions, the exponent 9 has been found to be 

B c = 0.13 ± 0.01 (5) 

and this value also has no measurable dependence on 
the initial (small) value of sublattice magnetisation m s (0) 
(Fig. 3). In agreement with the theoretical predictions 
H and in contrast with the simulation studies with global 



conservation algorithm M , this value of 9 has no numer- 
ically detectable dependence on the conserved magneti- 
sation too (Fig. 4) at least for too > 0.12. However, in 
agreement with both the studies gQ, the exponent has 
a different value for model A : 

9 A = 0.10 ± 0.01 (6) 

There is some indication that for very small values of Too 
(< 0.08), 9 might depend on mo (Fig. 4) but a detailed 
study in this region requires very large scale simulations. 
A few comments are in order. 

(i) Our estimates of 9 in 2d for model A is in agreement 
with previous estimates ||. The estimates for model A 
are obtained by putting too = in the present model (as 
too = corresponds to model A Jl2]|). 

(ii) That 9a — 9c in 2d indicates that the conservation 
is irrelevant here. This is because the specific heat ex- 
ponent a is negative here (|||, an d the estimate of 9c 
given by || is strictly true for 2 < d < 4 when the spin 
dimension n = 1. 

(iii) In 3d, our results 9a ^ 9c and independence of 9c 
of too are in qualitative agreement with the theoretical 
estimates ||. 

(iv) As regards the dependence of 9c on too in 3d, the 
discrepancy between our results and that of J7| is not sur- 
prising, because non-local conservation may be expected 
to give new values of exponents and can even change the 
universality classes. In fact, the numerical estimate of the 
dynamical exponent z in the globally conserved model C 
is found to correspond to that of model A, a result 
predicted by the analytical study of JnJ . 

(v) The results shown in the figures are for the largest 
sizes simulated. Simulation for smaller sizes, e.g., 199 x 
200 in 2d and 41 x 41 x 42 in 3d show that there is no 
detectable finite size effect. 

(vi) There could be some error in the estimate of 9 due 
to the error in the estimation of T c . We have, however, 
checked that 9 does not vary measurably when T c is var- 
ied about it's estimated value within the error bar. 

Lastly, we should mention that one can also estimate 
9 from the behavior of the autocorrelation function 

A(t,0) = (Si(t)Si{0)) - (Si(t)){Si(0)) ~ r x (7) 

where A = d/z — 6. In this method, one can put 
to s (0) = 0. One needs to know z very accurately to 
have an accurate estimate of 9 using this method. Also, 
a good estimate of 9 requires lattice sizes much greater 
than the ones simulated in the present study. Neverthe- 
less, we could calculate A(t, 0) and verify that A (and 
hence 9) is indeed independent of to . 

The authors are grateful to S. M. Bhattacharjee for 
valuable discussions and suggestions and to B. Zheng for 
useful comments. The computational facilities were pro- 
vided by CUCC. PS acknowledges financial support from 
DST project SP/S2/M-11/99. 
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Table 1 

Numerical estimate of critical temperature T c for 
different values of the conserved magnetisation mo- The 
estimates for 3d with toq ^ are in agreement with j7j. 



Dimension 


Size of Lattice 


m 




2 


1001 x 1002 





2.2692° 


0.1 


2.108 ±0.002 


0.2 


1.557 ±0.002 


3 


101 x 101 x 102 





4.5116" 


0.08 


4.453 ±0.004 


0.12 


4.376 ±0.004 


0.2 


4.142 ±0.004 


0.4 


3.009 ±0.004 



a Exact Result. b Approximate value, see e.g., [jjj 
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Fig. 1 : Early-time effect at critical temperature in two dimensions for a fixed value of the conserved magnetisation 
to , and various values of the initial sublatticc magnetisation to s (0). The lattice size was 349 x 350 and the results 
were averaged over 20,000 to 50,000 realisations of the system. 
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Fig. 2 : Early-time effect at critical temperature in two dimensions for a fixed value of the initial sublattice magneti- 
sation m s (0) and two values of the conserved magnetisation m . The line marked A is for Model A, i.e. to = 0. The 
lattice size was 349 x 350 and the results were averaged over 20,000 to 50,000 realisations of the system. 
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Fig. 3 : Early-time effect at critical temperature in three dimensions for a fixed value of the conserved magnetisation 
to , and various values of the initial sublattice magnetisation m s (0). The lattice size was 65 x 65 x 66 and the results 
were averaged over 5000 realisations of the system. 
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Fig. 4 : Early-time effect at critical temperature in three dimension for a fixed value of the initial sublattice 
magnetisation to s (0) and several values of the conserved magnetisation m . From top to bottom, the lines correspond 
to too =0.4, 0.2, 0.12, 0.08, (Model A); the results were averaged over 5000, 5000, 50000, 50000, 13000 realisations 
respectively. The lattice size was 65 x 65 x 66 in all the cases. 
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